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Abstract. We discuss two versions of the Frechet distance problem in 
weighted planar subdivisions. In the first one, the distance between two 
points is the weighted length of the line segment joining the points. In 
the second one, the distance between two points is the length of the 
shortest path between the points. In both cases, we give algorithms for 
finding a (1 + e)-factor approximation of the Frechet distance between 
two polygonal curves. We also consider the Frechet distance between 
two polygonal curves among polyhedral obstacles in TZ^ (1/ oo weighted 
region problem) and present a (1 + e)-factor approximation algorithm. 



1 Introduction 

Measuring similarity between curves is a fundamental problem that 
appears in various applications, including computer graphics and 
computer vision, pattern recognition, robotics, and structural biol- 
ogy. One common choice for measuring the similarity between curves 
is the Frechet distance, introduced by Frechet in 1906 [16]. The tra- 
ditional (continuous) Frechet distance 5f for two parametric curves 
P,Q: [0, 1] ^ 7^'^ is defined as 

5f{P,Q) = inf sup 5(P(a(r)),g(/3(r))), 

a,/3:[0,l]^[0,l] rg[o,l] 

where a and /3 range over all continuous non-decreasing functions 
with a(0) = /3(0) = and a{l) = /3(1) = 1, S* is a distance metric 
between points, and d > is the dimension of the problem. 

The Frechet distance is described intuitively by a man walking 
a dog on a leash. The man follows a curve (path), and the dog 
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Fig. 1. An example of the Prechet distance problem in weighted regions. 



follows the other. Both can control their speed independently, but 
backtracking is not allowed. The Frechet distance between the curves 
is the length of the shortest leash that is sufficient for the man and 
the dog to walk their paths from start to end. 

A reparameterization is a continuous, non-decreasing surjective 
function. The pair a, f3 of reparametrizations define how the end 
points of the leash, i.e. P{a{r)) and Q{(3{r)), sweep along their re- 
spective curves. We say that a, /3 is a matching between P and Q or 
a, /3 define a monotone walk from leash P{0)Q{0) to leash P{1)Q{1). 
In this We use s and t to denote a(r) and /3(r), respectively. We as- 
sume P and Q are polygonal, and denote by S{P) = {ui,U2, ■ ■ ■ , Up) 
and S{Q) = {vi,V2, ■ ■ ■ ,Vq) the sequences of vertices of P and Q, 
respectively. 

In this paper, we study the Frechet distance in weighted re- 
gions. In this problem, a 2D plane is divided into a subdivision 
R = {Ri, i?2) • • • ) Rn'} with a total of n vertices, with each region 
Ri & R having associated a positive integer weight Wi. The length 
of a path TT that stays within a region Ri of R is defined as Wi|7r|, 
where |7r| is the Euchdean length of tt. The length of a path in R is 
the sum of the lengths of the subpaths within each region of R. 

Let P{s)Q{t) be the leash with endpoints P{s) and Q{t). Let 
Ri{P{s)Q(t)) be the Euclidean length of the line segment with end- 
points P{s) and Q{t) within the region R^. We define the distance 
between two points P{s) and Q{t) on P and Q, respectively, as (a) 
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S{P{s),Q{t)) = Y:L,w,*R,{P{s)Q{t)) or (b) S{P{s),Q{t))= the 
length of the shortest (weighted) path from P{s) to Q{t). 

We formally define the Prechet distance problem in weighted re- 
gions as the following: Given a 2D weighted subdivision R and two 
parameterized polygonal chains P and Q in R, find the Frechet dis- 
tance between P and Q, where the distance between two points P{s) 
and Q{t) on P and Q, respectively, is defined either as in (a) or as 
in (b) above. 

Without loss of generality, we assume R is triangulated and P 
and Q lie on boundaries of R (see Fig. 1 as an illustration). We also 
discuss a special case in 7^^, where the weights are either 1 or oo 
(i.e., free space and obstacles). 

1.1 Motivation 

The motivation for studying such measurement comes from path/travel 
planning. Let us consider the following scenario. Suppose during a 
military operation, there are two teams of military units traveling 
on separate paths. However, before reaching their own destinations, 
the two teams want to maintain constant radio communications so 
that in case of any emergency one team can rescue the other team 
in time. The question is how should these two teams schedule their 
trips such that at any given time, the radio signal received by one 
team from the other is stronger than a threshold, or what is the 
optimal traveling schedule for each team such that minimum signal 
strength received is maximized. 

We use the Beer-Lambert law to model the decay (of the inten- 
sity) of an EM wave traveling through a region. Beer-Lambert law 
states that there is a logarithmic dependence between the transmis- 
sion T of the EM wave and the product of the absorption coefficient 
associated with the region a and the distance / traveled by the EM 
wave in the region, i.e. T = ^ — e~"', where and / are the ini- 
tial intensity (or power) of the EM wave and the intensity after the 
path, respectively [17]. If the EM wave travels through a set of re- 
gions R = {Ri, R2, . . . , Rk}, the transmission T can be written as 
T = eSi=i°*'% where «j and /j arc the absorption coefficient associ- 
ated with Ri and the distance traveled by the EM wave in R^, respec- 
tively. The decay of the wave intensity can also be expressed in terms 
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of the absorbance A which is defined as A = — log^^Q = J2i=i cah- 
Notice that if we treat i? as a weighted subdivision, where the weight 
of each region Ri is its absorption coefficient o;,, A is exactly the 
weighted length of the path traveled by the EM wave. Hence, this 
scheduling problem can be reduced to the Frechet distance problem 
in weighted regions. 

1.2 Previous Work 

The Frechet distance and its variants attracted considerable atten- 
tion in literature. Most previous work assumes an unweighted envi- 
ronment and can be divided into two categories, depending on the 
distance metric used. In the first category, the distance between two 
points is the Euclidean distance. In other words, the leash is always 
a line segment, as in case (a) above. Bending of the leash is not 
allowed. Frechet distances in this category are also referred to as 
non-geodesic Frechet distances. 

Although Frechet distance was introduced in 1906, it was not 
until 1993 that Alt and Godau [4] gave the first polynomial-time 
algorithm that computes the exact non-geodesic Frechet distance 
between two polygonal curves P and Q. They first gave an 0{pq) 
time algorithm for solving the decision version of the problem, which 
is the problem to decide, for a given constant Z\, whether the Frechet 
distance between two curves is at most A. Then, they showed how to 
solved optimization version of the problem, which is the problem to 
find the exact Frechet distance, in 0{pqlogpq) time by applying the 
decision algorithm and parametric searching technique. They also 
presented a more practical alternative solution, which takes 0{{p'^q+ 
q^p)^og{pq)) time, by running binary search over a set of critical 
values. 

Shortly after. Either and Mannila [15] introduced the concept of 
the discrete Frechet distance 6dF, which considers only order pre- 
serving pairing of vertices in the polygonal curves. They gave an 
0{pq) time algorithm for computing the discrete Frechet distance. 
They showed that the difference between the discrete Frechet dis- 
tance and exact Frechet distance between two curves is at most the 
length of the longest edge of the curves. As a result, discrete Frechet 
distance can be can be used for approximation of the exact Frechet 
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distance between two arbitrary curves. The approximation error de- 
creases as more vertices (Steiner pionts) are introduced to the two 
curves. 

Rote [20] explored the Frechet distance between more general 
curves. He gave an 0{pq logpg) time algorithms for finding the Frechet 
distance between piecewise smooth curves. The decision version of 
the problem can be answered in 0{pq) time. Wang et al. [7] studied 
the problem of measuring partial similarity between curves. They 
presented the first exact polynomial-time algorithm to compute a 
partial matching between two polygonal curves, which maximizes 
the total length of subcurves of them that are similar to each other, 
where similarity is measured in Frechet distance under Li or 
norm. Their algorithm takes 0{pq{p -\- qY\og{pq)) time. Buchin et 
al. [6] studied the Freechet distance between simple polygons and 
showed that it can be computed in polynomial time. They gave the 
first polynomial-time algorithm for computing the Frechet distance 
between simple polygons. 

In the second category, the distance between two points is the 
geodesic distance. Frechet distances in this category are also referred 
to as geodesic Frechet distances. The leash is allowed to bend to 
achieve to the minimum length. Maheshwari and Yi [18] study the 
case in which the curves lie on a convex polyhedron. By constructing 
a visibility diagram that encodes shortest path information for any 
pair of points on curves, one from each curve, they gave a Oi^ijP'q + 
pq'^)k'^ \og{pqk)) time algorithm, where k is the number of faces of the 
polyhedron. Cook and Wenk [14] gave a 0{k + log(A;M) logM) 
expected time or 0{k + log{kM)) worst-case time algorithm for 
computing the Frechet distance between two polygonal curves when 
the leash is constrained inside a polygon, where k complexity of the 
polygon and M is the larger of the complexities of the two curves. 
Chambers et al. [8] showed how to compute the homotopic Frechet 
distance between two polygonal curves in the presence of obstacles 
in polynomial time. The homotopic Frechet distance is the Frechet 
distance with an additional continuity requirement: the leash is not 
allowed to switch discontinuously, e.g. jump over obstacles is forbid- 
den unless the leash is long enough. 

Weighted region shortest path problems have been investigated 
in computational geometry for about two decades. Using Snell's re- 
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fraction low and the continuous Dijkstra algorithm, Mitchell and 
Papadimitriou [19] present an 0(n^log^^^) time algorithm, where 
N' is the largest integer coordinate of vertices and p is the ratio of 
the maximum weight to the mininiTim weight. Aleksandrov et al. [2, 
3] provide two logarithmic discretization schemes that place Steiner 
points along edges or bisectors of angles of regions in R, forming a ge- 
ometric progression. The placement of the Steiner points depends on 
an input parameter e > and the geometry of the subdivision. The 
(1 -|-e)- approximation algorithms in [2, 3] take -l-logn) log ^) 

and O(^logMog^) time, respectively. Sun and Reif [21] give an 
algorithm, called BUSHWHACK, which constructs a discrete graph 
G by placing Steiner points along edges of the subdivision. By ex- 
ploiting the geometric property of an optimal path, BUSHWHACK 
computes an approximate path more efficiently as it accesses only a 
subgraph of G. Combined with the logarithmic discretization scheme 
introduced in [2], BUSHWHACK takes 0(f (log \+\ogn) log i) time. 
Very recently, Aleksandrov et al. [1] gave a query algorithm that can 
find an e-approximate shortest path between any two points in 0{q) 
time, where ^ is a query time parameter. The preprocessing time of 
this algorithm is O ( ^^^J}j^ log ^ log^ i), where g is the genus of the 
discrete graph constructed by the discretization scheme. Cheng et 
al. [10] give an algorithm to approximate optimal paths in anisotropic 
regions, which is a generalized case of weighted regions. Their algo- 
rithm takes 0(^-^^n^log(^)) time, where p > 1 is a constant such 
that the convex distance function of any region contains a concen- 
tric Euclidean disk with radius 1/p. In weighted regions, the time 
complexity of the algorithm is improved to 0(^^^n^ log(^)) time, 
where p is the ratio of the maximum weight to the minimum weight. 
Very recently, Cheng et al. [11] also provided a query version of 
this algorithm that gives an approximate optimal patli from a fixed 
source (in an anisotropic subdivision) in O(log^) time. The prepro- 
cessing time is 0(^(log f 

1.3 Definitions and Preliminaries 

We denote by P{s)Q(t) the leash (line segment, also called link, 
or shortest path, depending of context) from P{s) to Q{t), where 

s,t e [0,1]. 
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Alt and Godau [4] introduced the free space diagram to solve the 
decision version of the non-geodesic Prechet distance problem (un- 
weighted case): Given two polygonal curves P and Q, and a positive 
constant A, determine if Sf{P, Q) < A. A free space diagram is a 
[0, 1] X [0, 1] parameter space such that each point (s, t) in the pa- 
rameter space corresponds to a leash with end points P{s) and Q{t). 
A free space cell C C [0, 1]^ is defined by two line segments, one 
from each curve. C corresponds to the leashes with endpoints on the 
two segments. The free space in the parameter space is defined as 
: S{P{s),Q{t)) < A}. That is, the free space represents all 
links shorter than A. A monotone path is a path monotone along 
both coordinate axes. There is a one-to-one correspondence between 
all possible matchings of P and Q and all monotone paths from (0, 0) 
to (1, 1) in the free space diagram. Hence, Sf{P, Q) < A if and only 
if there exists a monotone path in the free space from the bottom 
left corner to the top right corner of the parameter space. 

The discrete Frechet distance, introduced by Either and Man- 
nila [15], considers only the vertices of the (polygonal) curves. A 
coupling L between 6{P) and S{Q) is defined as 

such that ai = 1, bi = 1, am = p, bm = q, and for i = 1,2, ... ,m, 
we have Oj+i = or aj+i = ai + 1, and = 6j or = bi + 1. In 
other words, L is an order preserving pairing of vertices in P and Q. 
Note that each vertex can appear in L more than once. The discrete 
Prechet distance is defined as 

SdF{P,Q) = inf max S {uai ,Vb,). 

L i=l,2,...,m 

1.4 Our results 

Let e > be a positive constant given as part of the input. We 
also assume e < 1 at times. We have the following results: (1) 
For weighted regions in the plane, with S{P{s),Q{t)) — J27=i''^i * 
Ri{P{s)Q{t)) (case (a) above), we present a (l + e)-approximation al- 
gorithm that takes oipqN^log{pqN)) time, = 0(C(i?)(f (log i + 
log n) log i) is the total number of Steiner points used and C{R) is 
a constant associated with the geometry of the subdivision R; (2) 
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For weighted regions in the plane, with S{P{s), Q(t)) =the length 
of the shortest (weighted) path from P{s) to Q{t), we present a 
(1 + e)- approximation algorithm that takes 0(C(i?)^^(log^ + 

log J log^ ^) time, where ^ is a query time parameter related to 
computing shortest path and g is the genus of the graph constructed 
by the discretization scheme. (3) We give a (1 + e)-approximation 
algorithm for finding the Frechet distance between P and Q among 
polyhedral obstacles in TZ^ (1 and oo weights). To the best of our 
knowledge this is the first result for the TZ^ problem. The algorithm 
takes 0{C{Rfpq{l/e'^) log^(l/e)(n2A(n) log(n/e)/e^W log(n7) log(nlog7))) 
time, where 7 is the ratio of the length of the longest obstacle edge 
to the Euclidean distance between the two points, and A(n) is a very 
slowly-growing function related to the inverse of the Ackermann's 
function. 

2 The Line Segment Leash 

In this section, we discuss the case of the line segment leash. That 
is, S{P{s),Q{t)) = Y.l:LiWi*Ri{P{s)Q(t)). We first briefly discuss a 
pseudo polynomial exact algorithm for solving the decision version 
of this problem. Then, we address the optimization version of this 
problem directly and give a polynomial time approximation algo- 
rithm. 

2.1 An Exact Algorithm for the Decision Problem 

In this section, we give an exact algorithm extending Alt and Go- 
dau's algorithm [4] to solve the decision version of the problem. Re- 
call that the free space diagram introduced by Alt and Godau is a 
[0, 1] X [0, 1] parameter space such that each point (s, t) in the pa- 
rameter space corresponds to a leash with end points P{s) and Q{t). 
The free space diagram can be decomposed into 0{pq) free space 
cells such that each cell C corresponds to the leashes with endpoints 
on two segments, one from each curve, where p and q arc the number 
of vertices of P and Q, respectively. To determine if 6f{P,Q) < ^ 
for some positive constant A, we proceed as follows: 
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1. Partition each cell C into regions such that each region corre- 
sponds to leashes intersecting the same sequence of edges in R. 

2. For each region, find the free space with respect to the positive 
constant A by solving an 0{n) order bivariate polynomial system 
with 0{n) equations and inequalities. 

3. Determine if there exists a monotone path in the free space from 
(0,0) to (1,1). 

The partition in step 1 is an arrangement of 0{n) curves of the 
form st + cis + C2t + C3 = 0, where ci, C2, and C3 are constants. In 
the worst case, the total number of regions in the parameter space 
is (pqn^). Each region can be defined mathematically by 0{n) in- 
equalities. Once the free space in each region is computed, step 3 
can be solved by decomposing the parameter space and construct- 
ing a directed shortest path graph. We will address steps 1 and 3 in 
detail in later sections. Both steps can be solved in polynomial time. 
Next, we will show that Step 2 can be solved, which dominates the 
complexity of this algorithm. 

For each region, let f{s,t) be the weighted length of the leash 
P{s)Q{t) in terms of s and t, let and {ci, 62, ... , 6^} be the sequence 
of edges intersected by the leash. Note that k = 0{n) and Ci and 
Cfc are the two segments in P and Q, respectively. Following [9], 
the weighted length of leashes within the same region have the same 
functional expression, i.e. S{P{s),Q{t)) = Vl + J2tZi Wi \ - 

= vT+m^Ei-i , where ai,bi are constants, Wi is the 

m—mi I * ^i—L m—mi ' n i i i 
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weight of the region bounded by edges Cj and Cj+i, m, p are the 
slope and intercept of the leash, respectively, and rrii, Pi are the 
slope and intercept of e {ei, 62, ... , e^}, respectively. Let the end 

points of Ci be P(si) = {xi,yi) and P(s2) = (2^2,2/2), respectively, 
and the end points of Ck be Q(ti) = {xs,!/^) and Q{t2) — {x4,y4), 
respectively. For any s e [si, S2], we have 

P{s) = {xi+{x2-xi){s-si)/{s2-si),yi+{y2-yi){s-si)/{s2-si)), 
and similarly for any t & [ti, ^2], we have 

Q{t) = (x3 + (x4-X3)(t-tl)/(t2-il),y3+(y4-y3)(^-^l)/(i2-il)). 

Expressing m and p in terms of s and i , we obtain that 

CiS + C2t + C3 
TTl = 

dis + d2t + ds 

and 

c[st + C2S + c'^t + C4 



p 



dis + (i2^ + d^ 



where ci, C2, C3, c^, C2, Cg, c4, di, o?2, and o?3 are constants. It follows 
that 

f( ,\ _ /-I I / cis+C2t+C3 \2 v^fc ai{c[st+c'^s+c'^t+c'J/{dis+d2t+d3)+bi 
J{b,i) y ^ ^ \ dis+d2t+d3 ' ^j=l (cis+C2t+C3)/(<iis+ci2t+d3)-mi 

We can find the boundary of free space in this region by solving 
the following equation system: f{s,t) — A subjects to 0{n) inequal- 
ities that define the region. 

Since there are 0{n) fractional terms in f(s,t), we can convert 
this equation system into a bivariate polynomial system of degree 
0{n) which consists of one polynomial equation and 0{n) inequali- 
ties and requires pseudo polynomial time to solve [13]. 

We use Collins' cylindrical algebraic decomposition technique [13] 
analyzing the polynomial system, which takes 0{{fhn)^'^ d^) time, 
where m is the number of equation and inequalities in the polynomial 
system, n is the maximum degree of polynomials, r is the number 
of variables, A; is a constant and d is the maximum bit length of 
coefficients. Hence, step 2 takes 0{pqn^rfi'^^'>d'^^^^) time. 

This exact algorithm is not the one which ensures practical ap- 
plication. Next, we are going propose an approximation algorithm 
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for the optimization version of this problem. To approximate the 
Prechet distance, we discretize continuous space by placing Steiner 
points on edges of R. We place the Steiner points in such a way that 
all links can be grouped into sets. In each set, the difference between 
weighted lengths of any two links is small. Instead of finding the op- 
timal matching, we find the sequence of sets that the leash traversed 
in a specific matching. We then pick an arbitrary link from each set 
in the sequence, the maximum link length gives an approximation 
of the Prechet distance. 



2.2 Discretization Using Steiner Points 

We first discretize the continuous space by placing Steiner points in 
R extending the discretization scheme given in [2]. Let E be the set 
of all edges in R. Let V be the set of vertices in R. For any point 
V on an edge in E, let E{v) be the set of edges incident to v and 
let d{v) be the minimum distance between v and edges in E\ E{v). 
For each edge e E E, let d{e) — sup{d{v)\v e e} and let Ve be the 
point on e so that d{ve) = d{e). For each v & V, the vertex radius 
for V is defined as r(v) = — ^^-rr, where e is a positive real number 
defining the quality of the approximation, 5 is a lower bound on 
Sf{P,Q), and Wmax{v) is the maximum weight among all weighted 
regions incident to v. The disk of radius r{v) centered at v defines 
the vertex- vicinity oiv.B can be computed by setting weights of all 
regions to the minimum weight of R (assume all weights are greater 
than zero) and applying the continuous Frechet distance algorithm 
described in [4], where p and q are the number of vertices of the 
curves P and Q, respectively. 

For each edge e = V1V2 in E, we place Steiner points f j 1, f j 2, • • • , Vi^ki 
outside the vertex-vicinity of Vi, for i = 1,2, such that = r{vi), 

\vi,jVij+i\ = ed{vij), for j = 1, 2, . . . , fcj - 1, and Vi^k, = ^e- It fol- 
lows from [2] that the number of Steiner points placed on an edge is 
0(C(e)l/elogl/e), where C(e) = 0(j^log l^j ). Let TV de- 

note the total number of Steiner points and vertices of R. Since we 
set r(v) = — for each v e R, we have N = 0(C(R)(^(log i + 

logn) log i)), where C{R) = maxegij(^^^sgg^) is a function as- 
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Steiner edge 

ei / 



hourglass 




Fig. 3. An hourglass is the union of leashes intersecting with the same set of Steiner 
edges. 



sociated with the geometry of i?, and Wmaxi,^ is the maximum weight 
among all weighted regions incident to the end points of e. 

We refer to a line segment bounded by two consecutive Steiner 
points on an edge of /? as a Steiner edge. An hourglass is the union of 
all Icaslics (line segments) intersecting the same sequence of Steiner 
edges in the same order. Let H be an hourglass defined by a sequence 
of Steiner edges {ei, 62, . . . , 6^}, where ei e P and Ck & Q (See 
Fig. 3). 

Lemma 1. Let I and I' be two segments in H . Then, S{1) < (1 + 
2e)S{V) + 2eB. 

Proof. Let Ri^ , Ri^ , ■ ■ ■ , Rik-i weighted regions in if, such that 

Ri is between Cj and Cj+i. Recall that for a region R^ of R, Rk{l) 
denotes the Euclidean length of / within Rk. We have 

S{1) = J2w,^R,^{l) < + + 

j=l j=l 

where Wi^ is the positive weight associated with Ri^ . 

If a Steiner edge ej is outside of any vertex- vicinity, \ej\ < eRi._-^^{l') 
and \ej\ < eRi.{l'), due to the placement of the Steiner points. If 
Cj is incident to a vertex v, then \ej\ = r{v) = {tB) / {nWmaxiv)). 
One segment can intersect at most 0{n) Steiner edges inside vertex 
vicinities. The result follows. □ 
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Let Ih be an arbitrary segment in H, with endpoints on P and 
Q. Let a and p be two reparametrizations that define a matching be- 
tween P and Q. Let J — {Hi, H2, . . . , Hk} be the set of hourglasses 

traversed by the leash P{a{r))Q{P{r)). For an hourglass H & J, 
let Ih = {e\e = P{a{r))Q{/3{r)),r G [0,1], e G H}. Let H{a,/3) 
be the segment in Ih with the largest weighted length. That is, 
S{H{a,f3)) = maxee/^ S{e). Let 5{a,(3) = sup^gjo^^j S'(P(a(r)), Q(/3(r))). 

Lemma 2. | max/^g j S'(//f ) — S{a, (5) \ < 4e5(a, /3) . 

Proof. Applying Lemma 1, we have, 

S{Ih) < (1 + 2e)S{H{a, + 2eL', 

and 

S{H{a, /?)) < (1 + 2e)5(ij^) + 2eD. 
Assume e < 1/2, it follows. 



max5(/^^) > > (1 - 4e)<5(a, 

-Hfcj 1 -|- ze 



and 



maxSiln) < {1 + 2e) max S{H {a, p)) + 2eD < {1 + 4:e)5{a, p). 



□ 



We say that 5( J) = maxj^gj S{Ih) is a 4e-approximation of 5(q;, 
Given a sequence of hourglasses J = {-f^i, i^2, ■ ■ ■ , -f^fc}, we call J a 
legal sequence if there exists two reparmetrizations a, /3 that define 
a leash traversing the same sequence of hourglasses as J. 

Lemma 3. We can find a value 5'p{P,Q) such that \6'p{P,Q) — 
Sf{P,Q)\ < 4:edp{P,Q), that is, dp{P,Q) is a Ae- approximation of 
Sf{P,Q). 

Proof. Let a,j3 be the optimal reparametrizations that give the Frechet 
distance between P and Q. Let J be the sequence of hourglasses tra- 
versed by the leash defined by a and (3. Applying Lemma 2, we 
have, 

\S{J) - S{a, P) I = \S{J) - Sp{P, Q) I < 4e5(a, ^) = 4e5^(P, Q). 
So, there is a legal sequence of hourglasses that gives a 4e- approximation 

ofM^,g). □ 
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2.3 Prechet Distance Between Two Segments 

Next, for simplicity, we study a special case of the problem, in which 
each curve consists of one line segment. See Fig. 4 for an illustration. 
Similar to Alt and Godau's algorithm, we attack the problem in the 
parameter space D — [0, 1]^, where a leash P(s)Q{t) is associated 
with a point (s, t) e D. We refer to (s, t) & D as the dual point of 
the leash P{s)Q{t). 

Lemma 4. All leashes through a Steiner point v correspond to a 
curve in D, with equation Cy : st + CiS + C2t + c^ — 0, where C\, C2, 
and C3 are constants. 

Proof. Let v = (a, 6), P(0) = P(l) = (^2,^2), Q(0) = 

(2^3,2/3), = {xa.Va)- We have 

P(s) = {xi + (x2 - yi + (1/2 - 

Q{t) = + (0:4 - X3)t, y3 + (^4 - y3)t), 

and 

: 1/ - (z/i + {y2 - y.)s) = SftESEfeJfeE^lx - {x, + 

{X2 - Xi)s)). 

Since P{s)Q{t) passes through v, we obtain that 

{b-yi- {y2 - yi)s){xi -X3 + {x2 - xi)s - {x^ - X3)t)) = (a - xi - 

(x2 - xi)s){yi -y3+ {y2 - yi)s - {y^ - ^3)*), 

and thus we have 

Cy-.St^- CiS + C2t + C3 = 0, 

where Ci, C2, and C3 are constants. □ 

We call the curve the dual curve of Steiner point v. is 
continuous and monotone along both (s and t) axes. Obviously, v 
lies on a leash if and only if the dual point of the leash in D lies on 
Cy. Next we define the relative position of two leashes with respect 
to a Steiner point v. Given two leashes P{s)Q{t) and P{s')Q{t'), 
we say they are on the same side of v if and only if there exists 
two continuous functions a',/?' : [0, 1] — )■ [0, 1], such that a'(0) = s, 
a'{l) = s', /3'{0) = t, = t' and v ^ P{a'{r)Q{(]'{r)),yr G [0, 1]. 
Note that a', /3' are not necessarily monotone. Intuitively, two leashes 
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Fig. 4. A special case of the Prechet distance problem. 



are on the same side of f , if one leash can transform to the other 
one by sweeping its end points on their respective curves without 
crossing v. 

Lemma 5. divides D into two partitions, each partition corre- 
sponding to all links on the same side of point v. 

Proof. Given two hnks P{s)Q{t) and P{s')Q{t'), since there is a one- 
to-one correspondence between all possible walks from P{s)Q{t) to 
P{s')Q{t') and all paths from to {s',t') in D, if (s,t), 
are on different sides of C„, all paths between {s,t) and {s',t') must 
pass C„ at least once. That is, no walk exists between the two hnks 
without crossing v. □ 

We partition D by the dual curves of the Steiner points. Using 
the algorithm in [5], the partition can be computed in 0(A'"log N+k) 
time andO(iV + A;) space, where = 0(C(i?)(f (log i + logn) log^)) 
is the total number of Steiner points and k is the number of cells, 
which is 0{N'^) in the worst case. Let A denote the partition. 

Lemma 6. In A, each cell corresponds to an hourglass. 

Proof. If two points (s, t), (s', t') belong to the same cell, then there 
exists a walk between P{s)Q{t) and P{s')Q{t') that does not cross 
any Steiner point, i.e. P{s)Q(t) and P{s')Q(t') intersect with the 
same set of Steiner edges. If (s, t), (s', t') belong to different cells, 
then any walk between P{s)Q{t) and P{s')Q{t') must cross at least 
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one Steiner point. That is, P{s)Q{t), P{s')Q{t') do not intersect with 
the same set of Steiner edges. The result follows. □ 



The Frechet distance between P and Q can then be approximated 
as follows: 

1. Place Steiner points as described previously. 

2. Partition D by dual curves of all Steiner points. 

3. For each cell Z, choose an arbitrary leash I and assign S{1) as the 
weight of the cell. 

4. Find a monotone path T in D, from its left bottom corner, (0, 0), 
to its top right corner, (1, 1), such that the maximum weight of the 
cells traversed by T is minimized. 

Since the sequence of cells traversed by a monotone path in D 
corresponds to a legal sequence of hourglasses traversed by a leash 
following a match of P and and vice versa, by Lemma 3 the 
maximum weight of the cells traversed by T is a 4e-approximation 
of the Frechet distance between P and Q. 

2.4 Finding an optimal path in D 

We define the cost of a path between two points in D as the maximum 
weight of the regions traversed by the path. We decompose D by 
extending a horizontal as well as a vertical line from every vertex 
until it reaches the boundary of D. See Fig. 6 for an illustration. Let 
the new subdivision be D' . 

Lemma 7. Given an edge e in D' , let p be a point on e and let Tp 
be an arbitrary monotone path from (0, 0), i.e. the bottom left corner 
of D' , to p. Then, there exists a monotone path from (0, 0) to any 
point on e with the same cost ofTp. 

Proof. We consider e in D. Note that D consists of dual lines of 
Steiner points only. In we extend a vertical line as well as a hori- 
zontal line from the end points of e. Let Li, L2 be the two horizontal 
lines which define a horizontal slab and L^,^Li be the two vertical 
lines which define a vertical slab. Obviously, Tp has the same cost in 
D as it does in D' . Without loss of generality, assume enters the 
horizontal slab before it enters the vertical slab. Let v be the entry 
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point. Let S = {si, S2, ■ ■ ■ , Sk} be the set of curves on edges of D, 
which are bounded by the horizontal slab and traversed by Tp. Since 
the horizontal slab does not contain any vertex of D in its interior, 
no two curves in S intersect with each other. Hence we can construct 
a monotone path from v to any point on e that traverses S, i.e. we 
can contract a monotone path from (0, 0) to any point on e such that 
it has the same cost as Tp. See Fig. 5 for an illustration. The same 
argument holds if Tp enters the vertical slab first. □ 

□ 

Thus, given an optimal monotone path Tp from (0, 0) to an arbi- 
trary point p e e, we can construct a monotone path from (0, 0) to 
any point in e, which has the same cost as Tp and is also optimal. To 
find the optimal monotone path from (0,0) to (1,1), we construct 
a directed graph G having edges of D' as vertices. In G, a directed 
edge is added from a node Vi to a node V2 if, in D', Vi and V2 are 
on the boundary of the same cell and there exists a monotone path 
from vi to V2- The optimal monotone path can be found by running 
a modified Dijkstra's algorithm on G. 

Time complexity: The complexity of D' is N'^, where = 0(C(i?)(^ (log ^-1- 
logn) log ^)) is the number of Steiner points used, since each cell in 
D' has 0(1) edges and vertices, and thus each cell contributes 0(1) 
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edges to G. Then, G has 0{N*) vertices and 0{N*) edges and Di- 
jkstra's algorithm takes 0{N'^\ogN) time. 

2.5 Frechet Distance Between Two Polygonal curves 

We extend the algorithm above to approximate the Prechet distance 
between two polygonal chains P and Q. Recall that p and q are 
the number of vertices of P and Q, respectively. The parameter 
space D can be divided into {p — l){q — 1) subspaces, such that 
each subspace corresponds to all leashes bounded by the same two 
segments, one from each chain. Each subspace can be partitioned 
independently by introducing dual curves of Steiner points. The 
Frechet distance 5f{P,Q) can be approximated by finding an op- 
timal monotone path from the bottom left corner to the top right 
corner of D. The complexity of D is 0{pqN'^) and the complex- 
ity of the decomposed parameter space D' is 0{p^q^N^), where 
N = 0(C(i?)(^(log \ + logn) log \)) is the total number of Steiner 
points. It takes 0{p'^q^N^\og{pqN)) time to approximate 5f{P,Q)- 

Theorem 1. An approximation of the Frechet distance between two 

polygonal curves in a weighted subdivision can be computed in 
0{p^q^N^\og{pqN)) time, where N = 0(C(i?)(f (log i+logn) log i)) 
is the total number of Steiner points. 

3 Geodesic Frechet Distance 

In this section we study two versions of the geodesic Frechet distance 
problem. We do not require the leash to be homotopic, i.e. the leash 
is allowed to sweep discontinuously without penalty. For example, 
the leash can pass through or jump over obstacles. 

3.1 Geodesic Frechet Distance in Weighted Regions in 
TV 

Recall that the cost (weighted length) of a path in R is defined as 

the weighted sum of its Euclidean lengths within each region of R. A 
geodesic path between two points in i? is a path between those points 
that has minimum cost. Here, the distance between two points in R 
is the cost of the geodesic path between those points. 
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Fig. 6. Decomposition of D. 



We prove that the geodesic Frcchet distance can be approximated 
by the discrete geodesic Frechet distance. Given two polygonal curves 
P, (5, to approximate 5i?(P, Q), we need to add additional vertices to 
P and Q. We follow a similar approach as in [2], except that we define 
the vertex radius of a vertex f in P or Q as riv) = — ^^rr, where e is 
a positive real number defining the quality of the approximation, B 
is a lower bound on 6f{P, Q), and Wmaxi'v) is the maximum weight of 
regions incident to v. Recall that the disk of radius r{v) centered at v 
defines the vertex- vicinity of v. New vertices are placed on edges of P 
nd Q forming a geometric progression with ratios depending on e and 
on the geometry of R (see Section 3.1 for details). The total number 
of additional vertices introduced on each edge is 0(C(P)^ log^ 
where C{R) is a constant associated with geometry of R. Let P' 
and Q' be the new polygonal chains and let p' — 0{C{R)^log^ ^) 
and q' = 0{C{R)^log^ ^) be the number of vertices of P' and Q', 
respectively. 

Lemma 8. The discrete geodesic Frechet distance between P' and 
Q' gives an e- approximation of the geodesic Frechet distance between 
P and Q, I.e. (1 - e)5i.(P, Q) < 6dF{P': Q') < (1 + e)5F(P, Q). 

Proof. Here, d{a^ b) denotes the weighted length of the line segment 
ab. Let a, f3 be two optimal reparametrizations between P' and Q'. 
Obviously, they are also an optimal matching between P and Q. Fol- 
lowing [15], we define a coupling L between S{P') = {ui, U2, . . . , Up/} 
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and 5{Q') = {vi,V2, ■ ■ ■ ,Vq'}. For each point u G S{P'), let t{u) be 
the smallest value such that P'{a{t{u))) = u, and for each point v 
in 5{Q'), let s{v) be the smallest value such that Q'{^{s{v))) — v. 
First, we add {ui^vi) to L. Let i and j be the number of distinct 
points in 5{P') and 5{Q')^ respectively, added to R. If j — q' or 
< 5(^^+1), then add {uiJ^i^Vj) to L. We have 

d{ui+^,Vj) < S(P'(a{t(ui+i)))Q'(/3(t(ui+i)))) + d(vj,Vj+i) 

and 

S(P'(a(t(ui+i)))Q'(/3(t(ui+i)))) < d(ui+i,Vj) + d(vj,Vj+i). 
If VjVj+i is not contained in any vertex vicinity, we have 

d(vj,vj+,) < eS(P'(a(t{ui+i))),Q'(p{t(ui+^)))) < e5F(P,Q). 
If VjVj+i is inside a vertex vicinity, we have 

d{vj,Vj+i) <eB< €6f{P,Q). 

We obtain that 

\S{P'{a{t{ui+^))),Q'{m^i+lm - d{u,+,,vj)\ < eSpiP.Q). 

li i = p' or t{ui+i) > s{vj+i), then we add {ui, vj+i) to L. Similarly, 
we have 

\S{P'{a{s{vJ+^))),Q'{p{s{vJ+^m - d{ui,Vj+,)\ < e5F{P,Q), 
and thus 

(1 - e)5FiP, Q) < 5dF{P', Q') < (1 + e)5p{P, Q) 

□ 

Instead of finding an exact geodesic path between two points, 
we can apply existing shortest path algorithms in weighted regions, 
which find an approximate path. Let 6'^p{P', Q') be the approximate 
discrete Frechet distance computed by replacing the exact shortest 
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path algorithm by the approximation algorithm, which gives an e- 
approximation of the geodesic distance between points. We have 

S'dFiP',Q') < {l+e)SdF{P,Q) < (l + e)%(P,Q) < (l + 3e)5p(P, Q). 

Similarly, assuming that e < 1/3, we have 

S',^{P',Q')>{l-3e)6F{P,Q). 

Thus, we obtain 

\SdFiP',Q')-MP,Q)\ < Se6p{P,Q). 

If a non-query based approximation algorithm is used, 6'^p{P',Q') 
can be computed in time 0{p'q'T{n, e)) = 0{C{Rf f{\og^ i)r(n, e)), 
where T{n, e) is the time to compute an approximate shortest path 
between two points. If a query-based approximation algorithm is 

used, then 6'^p{P' ,Q') can be computed in 

0{p'q'QUERY{n,e)+PRE{n,e)) = 0{C{Rf^{\og^ l)QUERY{n, e)+ 
PRE{n,e)), where QUERY{n,e) is the query time and PRE{n,e) 
is the preprocessing time of the algorithm. For example, if we use the 
algorithm proposed by Aleksandrov et. al. [1], the algorithm takes 
0(C(i?)2f(log^i)g+ (^^logMog^i) time, where g is a query 
time parameter and g is the genus of the graph constructed by the 
discretization scheme. 

Theorem 2. A 3e- approximation of the Frechet distance between 
two polygonal curves in a 2D weighted subdivision R can be com- 
puted in 0(C(i?)^P^(log^ 7)^+ ^^^2}})^ log 7 log*^ ^) time, where q is a 
query time parameter, g is the genus of the graph constructed by the 
discretization scheme, C{R) is a constant associated with geometry 
of R, and n is the number of vertices of R. 

3.2 Geodesic Frechet Distance in TZ^ with Obstacles 

In this subsection, we briefly discuss the geodesic Frechet distance 
problem in 1 or cxd weighted regions in (that is. among obstacles 
in TZ^). Let i? be a weighted subdivision in 7^"^ with a total of n 
vertices. The weight of each region Ri E R is either 1 or cxo. Given 
two polygonal curves P and Q in R, we want to find the Frechet 
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distance between P and Q, where the distance between two points 
in R is defined as the length of the geodesic path between those 
points, i.e. the length of the shortest obstacle-avoiding path. 

Wc set r{v) = eB and add additional vertices on P and Q as 
described in Section 4.1. Let P' and Q' be the new curves. 

Lemma 9. The discrete Frechet distance between P' and Q' gives 
an e- approximation of the Frechet distance between P and Q, i.e. 
{l-e)SF{P,Q)<SdF{P',Q') < il + e)5F{P,Q). 
Proof. Similar to the proof of Lemma 8. □ 

Let 6'^p{P',Q') be the discrete Frechet distance computed by a 
shortest path approximation algorithm, which gives an e- approximation 
of the shortest path. We have 

\6',AP',Q') - 6f{P,Q)\ < 3e6F{P,Q). 

S'apiP', Q') can be computed in 0{C{Rfpq{l/€^) log^(l/e)T(n, e) time, 
where C{R} is a constant depending on the geometry of the problem 
and T{n, e) is the time to approximate the shortest path between 
two points in R. For example, we can use the approximation algo- 
rithm given by Clarkson [12], which takes 0(n^A(n) log(n/e)/e'^ + 
log(n7) log(nlog7)) time, where 7 is the ratio of the length of 
the longest obstacle edge to the Euchdean distance between the 
two points, and X{n) is a very slowly-growing function related to 
the inverse of the Ackermann's function. Thus, our algorithm takes 
0{C{Rfpq{l/e^) log^(l/e) {n^X{n) log(n/e) /e^+n^ log{n-f) log(n log 7))) 
time. 

Theorem 3. A ?>e- approximation of the geodesic Frechet distance in 
TZ^ with Obstacles can be computed inO{C{RYpq{l/ 1^) log'^(l/e)(n^A(n) 
log(n/e)/e^-|-n^ log(n7) log(nlog7))) time, where^ is the ratio of the 
length of the longest obstacle edge to the Euclidean distance between 
the two points, and \{n) is a very slowly- growing function related to 
the inverse of the Ackermann's function. 

4 Conclusion 

In this paper, we discussed three versions of the Frechet distance 
problem in weighted regions and presented an approximation algo- 
rithm for each version. First, we discussed the non-geodesic Frechet 



22 



distance problem in planar weighted regions. We showed that we 
can approximate the Prechet distance by using a parameter space, 
D, where each leash is associated with a point in D, and constructing 
a discrete graph G from D. We then discussed two geodesic Frechet 
distance problems, in planar weighted regions and in 1 or oo weighted 
regions in TZ^, and showed that in both cases, by adding additional 
vertices to the polygonal curves, the discrete Prechet distance can 
be used to approximate the continuous Prechet distance. 
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